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ABSTRACT 

We present Monte Carlo studies of charge expectation values and charge fluctuations 
for quasi-particles in the quantum Hall system. We have studied the Laughlin wave func- 
tions for quasi-hole and quasi-electron, and also Jain's definition of the quasi-electron wave 
function. The considered systems consist of from 50 to 200 electrons, and the filling frac- 
tion is 1/3. For all quasi-particles our calculations reproduce well the expected values of 
charge; —1/3 times the electron charge for the quasi-hole, and 1/3 for the quasi-electron. 
Regarding fluctuations in the charge, our results for the quasi-hole and Jain quasi-electron 
are consistent with the expected value zero in the bulk of the system, but for the Laughlin 
quasi-electron we find small, but significant, deviations from zero throughout the whole 
electron droplet. We also present Berry phase calculations of charge and statistics pa- 
rameter for the Jain quasi-electron, calculations which supplement earlier studies for the 
Laughlin quasi-particles. We find that the statistics parameter is more well behaved for 
the Jain quasi-electron than it is for the Laughlin quasi-electron. 
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1 Background, motivation and summary of main re- 
sults 



The perhaps most striking feature of the fractional quantum Hall effect is the existence of 
quasi-particles with fractional charge and statistics. There now exists direct experimental 
evidence for the existence of fractionally charged particles , but much of the evidence for 
fractional statistics in the quantum Hall system is still indirect and based on the theoretical 
description of the effect. 

For the plateau states corresponding to filling fractions i/ = 1/m, m odd, the many- 
electron wave functions introduced by Laughlin are generally accepted as giving an 
essentially correct representation of the true quantum states (under the idealized condition 
of a homogeneous background potential). Also the quasi-hole excitations are described by 
simple wave functions, and there are robust arguments for these to be essentially correct. 
The charge and statistics of the quasi-holes have been determined by Berry-phase calcula- 
tions and agree with the claim that they are anyons |0. The situation is less clear for the 
quasi-electrons. Different trial wave functions have been suggested, notably by Laughlin 
and Jain but the conclusions based on these are less convincing than those for the 
quasi-holes. 

Although the indirect evidence for the charged excitations of the fractional quantum 
Hall effect to be anyons is rather convincing, a more direct verification would certainly 
be interesting. It is a challenge to establish experimental evidence for the (fractional) 
statistics of the charge carrying excitations in a similar way as their fractional charge has 
been found. However, also on the purely theoretical side it is of interest to examine further 
the anyon aspects of the quantum Hall effect. 

In a previous paper ^ we examined the anyon representation of the Laughlin quasi- 
holes in some detail. In that paper also the quasi-electrons were discussed, although much 
more briefly. In particular it was pointed out that in order to reproduce expected values 
for charge and statistics parameters by Berry phase calculations, approximations had to 
be done which could not readily be justified. This motivated subsequent numerical studies 
of quasi-electrons, by Kj0nsberg and Myrheim [Q, for systems with up to 200 electrons. 
For the charge parameter bulk values were found which were close to the expected value. 
However, for the statistics parameter the rather surprising result was that no stable bulk 
value was found. The numerical studies were done by applying the quasi-electron wave 
functions introduced by Laughlin. These results were quite different from the corresponding 
numerical results obtained for quasi-holes. 

The results obtained in Ref. have motivated the present work, where we examine 
charge expectation values as well as charge fluctuations for quasi-holes and quasi-electrons.Q 
The fractional charge of the physical quasi-particles is expected to be a sharply defined 
quantum number, and it is of interest to check whether this is true for the suggested 

^We are obliged to Hans Hansson for his suggestion to follow up the results of the Berry phase calcu- 
lations by examining the charge fluctuations. 
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trial wave functions. There exist some general arguments which link fractional statistics 
to fractional charge 0, and although the precise conditions for this to be true are not 
clear, they suggest that the problems which are seen in Berry phase calculation of the 
statistics parameter may also show up in the charge calculations. We have in particular 
been interested in examining the possibility of long range fluctuations in the quasi-electron 
charge. 

The calculations of expectation values and fluctuations of charge have been done by 
the Monte Carlo method for quasi-holes and quasi-electrons corresponding to the m = 3 
state. We have in particular compared results obtained with Laughlin's and Jain's 
definitions of the quasi-electron wave functions. The systems considered have electron 
numbers varying from 50 to 200. The main results of the calculations are the following: 

For quasi-holes the expected bulk values of the charge and the charge fluctuations 
are well reproduced. The expectation value of the charge is —1/3 (times the electron 
charge) and the charge fluctuation is zero. Numerically small statistical fluctuations are 
present, but no significant deviations from these values. For Laughlin's quasi-electron 
the numerical results for the charge expectation values are consistent with the expected 
value 1/3. However, for the charge fluctuations we find small, but statistically significant 
deviations from the value zero in the bulk. This is the case even for an electron number of 
200. For Jain's definition of the quasi-electron wave function we again reproduce expected 
results, 1/3 for the charge and vanishing fluctuations. Differences between these values 
and the numerical results are within small statistical errors. 

These results indicate that the problems earlier found in Berry phase calculations of the 
statistics parameter may be an artifact of Laughlin's special definition and not a signal 
of long range effects for a generic quasi-electron state. To check this more directly we have 
also examined the charge and statistics parameters found from Berry phase calculations 
of the Jain quasi-electron state. We find for this wave function a stable bulk value of the 
statistics parameter consistent with the value —1/3. 

Throughout the paper, we are using dimensionless complex coordinates z = [x + 

iy), with £b = 1/V eB, and we set h, c = 1, B being the magnetic field. 

2 Quasi-particle charge and charge fluctuations 

In general it is a subtle problem to define localized and sharp charges in quantum many 
body systems. The naive definition of a charge operator which gives the charge within a 
finite region A of space is (in 2 dimensions) 

Q{R)= f d^xp(f,f). (1) 

J A 

where p{x,y) = ip\x)ip{y) is the single particle density operator and A is taken to be a 
circular area of radius R. The charge measured relative to the ground state is defined by 
a trivial subtraction 

C{R) = Q{R)-m{Rm , (2) 
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where |0) is the many body ground state. When the ground state is represented as a filled 
Fermi sea, the relative charge C{R) can alternatively be defined by normal ordering the 
density operator with respect to the ground state. 

States with a well-defined particle number are eigenstates of the total subtracted charge 
operator C{oo)] in particular C*(oo)|0) = 0. For the charge operator corresponding to a 
finite area A, this is not the case, even if the radius R is taken to be very large. This can 
be seen by considering the charge fiuctuationQ 

{AQ{R)f = (QiRf) - {Q{R))' , (3) 

which does not vanish even for the ground state. In a relativistic field theory the charge 
fiuctuation in fact diverges due to contributions from particle-antiparticle pairs of arbitrary 
high momenta. In a non-relativistic theory the fiuctuation is finite, due to the finite depth 
of the Fermi sea, although it will in general be large. The case we consider here is even 
more well behaved, since the particles are restricted to the lowest Landau level. A deep 
Fermi sea would correspond to a situation with many filled Landau levels. 

For states with short range correlations, we expect the fiuctuation to be an edge effect, 
and to demonstrate this more explicitly we rewrite the charge fiuctuation in the following 
form, 

{AQ{R)f = - f d?x f d'yc{x,y) (4) 

where A'^ , is the complement of the area A, and c(x, y) is the density-density correlation 
function defined by, 

c{x,y) = {p{x)p{y)) - {p{x)){p{y)) . (5) 

For a homogeneous ground state the correlation function only depends on the relative 
distance, c{\x — y\) = c(r), and with an exponential fall off for large r the integrals in 
Eq. (Q) only get contributions close to the boundary of A and A'-' . This gives for the 
charge fluctuation a linear dependence on the radius R, when R is much larger than the 
correlation length. 

For the incompressible states of the quantum Hall system, the density correlations are 
short range, with a correlation length of the order of the magnetic length is- To illustrate 
the form of the charge fluctuation {AQ{R)y in this case, we consider the case of a flUed 
lowest Landau level, and from now on all lengths are measured in units of ^/2iB and 
integrals are given on dimensionless form. In the limit of an inflnite system we then have 
an analytic expression for the correlation function, 

c(r) = -(5(f)--e-'-') , (6) 

TT V IT J 

and the charge fluctuation can easily be calculated. This is the case even for a flnite 
system, and in Fig. |I| the functional form of (A(5(-R))^ is shown for an electron droplet 

^ Clearly, the subtraction in (Q) does not change the fluctuations, i.e. (AC(i?))^ = (AQ(i?))^. 
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corresponding to 50 electrons. The linear dependence is clearly seen for an intermediate 
range, where R is larger than ^B but smaller than the size of the system. For the Laughlin 
states we expect a similar behaviour. In this case we can use Laughlin's plasma analogy to 
demonstrate the exponential fall-off of the correlation function, a behaviour corresponding 
to the screening of charges in the plasma. Earlier numerical calculations of the correlation 



function |]T2[ also show the exponential fall-off, and we have by direct numerical calculations 
found a functional form of (AQ(i?))^ for the z/ = 1/3 state which is similar to the one shown 
in Fig. 0. 




Figure 1: The charge fluctuation (AQ(i?))'^ for a system with 50 electrons and filling factor m = 1. The 
two curves show the fluctuations of the ground state and of the state for one quasi-hole. 



It is a non-trivial problem to construct a local charge operator which is insensitive 
to the background fluctuations of the ground state. In the context of a one-dimensional 
fermi system with fractionally charged solitons, the problem was solved by Kivelson and 
Schrieffer P] and Bell and Rajaraman who defined a charge operator using a smooth 
spatial cutoff. However, in 2 (and higher) dimensions this is not sufficient, also a cutoff 
in energy is needed, as discussed by Goldhaber and Kivelson |jlO|. In principle such a 



definition could be used also for the fractional quantum Hall states, however to implement 
this in the form of an explicitly defined charge operator seems difficult. 

In this paper we have not made any attempt to introduce such a redefined charge 
operator. Instead we have used, as the relevant measure of the sharpness of a localized 
charge, the subtracted quantity 



D{R) = AQ{R) - AQo{R) 



(7) 



where AQq{R) is the fiuctuation of the ground state. Each of the two terms in this 
expression are finite, and the difference is therefore well defined. The use of this expression 
is based on the intuitive idea that if an excited state can be characterized by a sharply 
defined, localized charge, that means that the charge fiuctuation rapidly approaches the 
ground state value when the integration area A is extended outside the region which 
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characterizes the size of the excitation. In fact, it seems reasonable to assume that if the 
fluctuation vanishes for any reasonably defined local charge operator, then also the simple 
difference (|^) will vanish. The only problem with the approach used here is of numerical 
character. The fluctuation will be a small number calculated as the difference between 
large numbers, and this put strong limits on the accuracy of the numerical calculations. 

For the sake of illustration we show in Fig. |^ the charge expectation value and the 
charge fluctuation, with such a subtraction included, for the special case of a Laughlin 
quasi- hole at the integer filling factor 1. In this case all the quantities needed in order to 
calculate the functions C{R) and D{R) can be found analytically^, and we show the results 
for a system of 50 electrons. The excess charge, here negative, exponentially builds up to 
the plateau value (—1) as R increases from zero, and the charge fluctuations vanish once 
the plateau is reached. We also see edge effects at the droplet boundary R ^ 7.1. The 
figures clearly show that our definitions in Eqs. (H) and (|^) reveal the basic properties of 
interest. 




Figure 2: Quasi- hole charge Ch(-R) and charge fluctuation D\^{R) for a system with 50 electrons and 
filling fraction 1/m = 1. 

Before concluding this section on general definitions and relations, we shall introduce 
a quantity that will be very useful in the coming calculations. Consider 

F{R) = f d\id\p2{u,v), (8) 

J A 

where p2{u,v) is the diagonal two-electron density operator. Using 

p{u)p{v) = p2{u, v)+5{u- v)p{u, v) (9) 

■^The wave functions will be given in full detail in section ^, but for now our point is only that of 
illustration. 
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we can then express the fluctuation of the operator Q as 



{/\Q{R)f = {Q\R)) - {Q{R))' (10) 
= Q{R) + F{R)-{Q{R)f. (11) 

The functions Q, F are expectation values of Q,F in the specific state. In terms of the 
iV-electron normalized wave function for the state they are 

Q{R) = N j^d\i j (fz2---(fzN I ipiu,Z2,---,ZN) 1^ (12) 
F{R) = N{N - 1) / d\Ld% [ d% ■ ■ ■ d^N \ ^(m, v.z^,---, zn) ? (13) 



The expression (0) will be used when investigating the charge fluctuations of the Laughlin 
quasi-electrons, whereas ([T0|) will be used in the case of Laughlin's ground state and quasi- 
hole state as well as for Jain's quasi-electrons. 

3 Numerical methods to compute Q and AQ 
3.1 General introduction 

In all our computations we have used Monte Carlo integration with importance sampling 



according to the Metropolis algorithm |Tl|, O], and a brief discussion of some points of 



special importance for our calculations are given in appendix |^. The specific probability 
distribution used to generate electron configurations varied from case to case. Each of them 
will be given in the sections following below, where we give descriptions of our methods 
that are detailed enough to allow the calculations to be reproduced. 

Notice that in order to find the quasi-particle charge fluctuation D{R) we need to find 
the two numbers AQ{R) and AQq{R), and then the difference between them. In general, 
this difference is very small compared to the numbers themselves. Since our aim is to 
decide if the difference is significantly different from zero, both AQ{R) and AQo{R) must 
be known with great precision. This is itself a challenging problem since each of them in 
turn is found from a difference between the numbers {Q'^{R)) and {Q{R))'^ which are large 
for large R. For all examined cases except for Laughlin's quasi-electrons, this problem 
has been solved by evaluating (Q^) and (Q)^ simultaneously, and benefit from the fact 
that the statistical errors then are correlated and tend to cancel. For the Laughlin quasi- 
electrons this method is inappropriate. Instead we have found values for Q and F defined 
in Eqs. (|12D and ([l^) . In the following discussion we will use the subscripts "h, Le" and 
"Je" to refer to quasi-holes, Laughlin quasi-electron and Jain quasi-electrons, respectively. 
The subscript "0" refers to the ground state. 
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3.2 The Laughlin ground state and quasi-hole state 

The definitions (|^) and (0) of tlie cliarge and cliarge fluctuations for tlie quasi-particle, 
require knowledge of the mean value and fluctuations of the ground state charge. For 
electrons at filling fraction 1/m the ground state is given by Laughlin's wave function 

^o(;2i, ■■■z*j,) = e-^Sf.J^^P Y[(z^ - z^r, (14) 

k<i 

This describes an electron droplet with uniform particle density (in the bulk) and a radius 
approximately \fmN |jT2|, ITSf. In this case the numerical computations are simple. We use 
the true electron density to generate the electron configurations (2:10,, ■ ■ ■ , -Zatq), that is we 
generate events according to the probability density, 

Pq{zi,---z*^) = ^{ipoizi, ■ ■ ■ zIj)\^ , (15) 
-'0 

where 

' d'^'zlM^^-.-z^j,)]' (16) 



is the normalization integral. For each specific configuration a the charge Qa inside the area 
A = 7tR^ is found simply by counting the z:s that satisfy \z\ < R. From n configurations 
the expectation values of Q{R) and Q^{R) are estimated (simultaneously) by 

in 1 

(o|Q(i?)|o) = - E m'{R)\o) = -EQl (17) 

and the ground state fluctuation AQq{R) is then found using Eq. (plQl). 

The same method can also be used for a Laughlin quasi-hole, which is given by the 
wave function 

k<l 1=1 

This expression describes an electron droplet with an electron density which is the same as 
in the ground state, except for a small depleted region around the point 2:05 the position of 
the quasi-hole [|12|. Again we choose the probability distribution used to generate electron 
configurations to be the (normalized) electron probability density itself; 

Py,{zi, ■ ■ ■ z]^) = ^\ipY,{zi, ■ ■ ■ z*j^)\'^ , where h= ( d^^z\ipY,{zi,- ■ ■ z]^)\^. (19) 

For Zo = the expectation values (h|Q(-R)|h) and (h|Q^(-R)|h) are then estimated, using 
the procedure described above. 
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3.3 The Laughlin quasi- electron state 



The simple numerical method described above can not be used in the case of a Laughlin 
quasi-electron. The reason for this can immediately be seen from the wave function, which 
reads 

^Le(^i, ■■■z%)= e-^S.-l^^l n(5.. - ^o) Ui^k - zir. (20) 

i=l k<l 

In this case the wave function contains derivatives, and these must be evaluated analyt- 
ically before the expressions can be used in a computer program. Although in principle 
straightforward, we know of no efficient way of doing this for sufficiently many electrons^. 

To avoid this problem, we have considered two different methods for computing the 
charge and charge fluctuations of a Laughlin quasi-electron. The first uses "brute force", 
and the numerical convergence is slow, although much better for small than for large R. 
The second method converges faster, and is best fitted for the bulk. Hence, the domains of 
validity of the two methods are complementary, which allows us to calculate for the whole 
range of R. Also there is some overlap where we can check against one another the results 
obtained by use of the two methods. 



Method 1 

We will here present the method most appropriate for small R. Consider then the expec- 
tation value of the charge operator when the quasi-electron is put at the origin; 

QUR) = (Le|Q(i?)|Le) (21) 

= f / d^.i/d^(^-^).|V^Le(.i,---,.^)r (22) 

= f / d\ [ d^(^-).e-i:tj-i^n I ^^-^^ r m\ p -i) i r(23) 

iLe J A J j^^i i=2 ^1 ^« 

where /lb is the normalization integral, and is the integrand in Eq. (^31). The expression 
in Eq. (^) is obtained by partially integrating the coordinates Z2, - ■ ■ , z^, but differenti- 
ating the VanderMonde determinant directly with respect to Zi, which is integrated over 
the finite area A. We emphasize that p'l^^ is not the true electron density of the system, in 

That the Laughlin quasi-electron represents challenging numerical problems is well known. Several 
authors have calculated the single electron density with a quasi-electron placed at the origin of a circular 
droplet, and with results that are not in total agreement with one another. The discrepancies are restricted 
to the behaviour close to the origin. Haldane and Rezayi and later Morf and Halperin found, as 
opposed to Laughlin jl^, that the single electron density has a dip at the origin. In Ref. jl^ this dip even 
drops below 1. 
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fact, since it takes on negative as well as positive value, it cannot be taken as a probability 
density. Nevertheless, the normahzation integral Ilg can be expressed in terms of p'^^^, 



,Le= / d2^^|^Le(^l,---,^^)|' = J d'^'z p'^^{z^, ■ ■ ■ , Z%) , (25) 

and this enables us to write, 

^ . . _ /^d^Zi/d^(^"%pL(^,^*) _ ,, /Ad^^i/d^^^-^^^bLe(^,^*)|sgn(pL) ...^ 

^^^^ > Jd^^ZpU^,Z*) /d2%bL(^,^*)|sgn(p'Le) ' ^ ' 

where sgn(pLe) is the sign function. Eq. (l26| ) tells us that if the electron coordinates are 
generated according to the absolute value Ip'i^^^z, z*) \ then the expectation value QheiR) can 
be estimated as the ratio between the Monte Carlo expectation values of S'(-R— l^i |) sgn(pLe) 
and sgn(pLg), where S{x) is 1 when x > and otherwise. In practice we thus do as 
follows: For each electron configuration a determine whether particle 1 is inside R. If it 
is, set ta = ±1 according to the sign of p'^^{z,z*). Otherwise set = 0. Independent of 
the position of zi set Ua = ±1 according to the sign of p'-^^. The expectation value is then 
found by 

QUR) = ^F^. (27) 

Notice that the function |pLg(-2,z*)| is not symmetric in all variables since it treats Zi in 
a special way. This implies that the numerical convergence rate is times slower than it 
would have been if we could use a symmetric function, everything else being equal. 

Since the electron coordinates were not generated according to the true electron density, 
the expectation value {Q^{R)) can not be found using the simple method described for the 
Laughlin ground state and quasi-hole. Instead we turn to Eq. (pA| ) and compute 

FLe(i?) = (Le|F(i?)|Le) (28) 

= " / d\d% [ d'^^-% |^Le(^l, ■ ■ ■ , ^^)r 



(29) 

^^^^ / d\d\ f d'(^-'he-^i.\^^\' n \^>^ - ^irmi^ki' - 1) (so) 



n 

k<l k=3 

N N 

X 



i^l ^1 - ^2 - Zj [Zi - Z2) 

The VanderMonde determinant has now been differentiated with respect to both zi and 
Z2, while the other coordinates were integrated by parts. The integrand in Eq. ( ^I]) is more 
complicated than in Eq. (|2^). However, if we now define integrand of 

Eq. (0), then the normalization integral Jlg can still be expressed as = / d'^^zp'^^{z, z*), 
and the value of Fi^e{R) be estimated as the ratio between two Monte Carlo estimates, 
analogously to the method above. Notice that in the present case the function \p'^^{z,z*)\ 
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treats two coordinates differently from the others. The convergence rate of this non- 
symmetric treatment will be A^(A^ — 1) times slower than a symmetric one, again everything 
else being equal. 

The slow convergences for large N that we have referred to, imply problems for the 
numerical calculations discussed here. Nevertheless, for small values of R we do obtain well 
converged results. These can be used as a complement to the bulk results found by use of 
the method of the next section. In addition, the ranges of validity of the two methods do 
have some overlap, and we hence check the results against one another. We also hope that 
the treatment of "negative probability densities" is instructive, and would like to mention 
that also Ref. |T^ contains a discussion of this topic. It will be used again in the next 
section. 



Method 2 

This section presents a method that improves Method 1 in two important ways: First, a 
function symmetric in all electron coordinates is used to generate configurations, and thus 
we achieve a huge improvement in the convergence rate. Second, the same electron config- 
urations are now used for calculating both Q^eiR) and FLe(-R). However, to determine the 
functions we need to use numerical derivatives in addition to the Monte Carlo integration. 
Define the quantity 

N ^ 

j=l k<l 

Of course, pLe is nothing else than the integrand appearing if we use integration by parts 
to all coordinates in Eqs. (12^, 123), and it is clear that 



'Le 



d'^'z |V^Le(^l, ■■■,Z*^)\'= d'''zpUzu ■ ■ ■ , Z*n)- (33) 



Q{R) = ^ f d'z, f d2(^-i)^pLe(^i, ■ ■ ■ , ^^v), (34) 



In addition, define the auxiliary functions, 

N r 

he 'a 

F{R,,R,) = ^^^~^^ [ d'zJ d'zJd'^''-'hpUzu---,z*N)- (35) 

JLc "'^1 "'-42 

Notice that the areas Ai and A2 in general are different. 

The auxiliary functions Q{R) and F{Ri, R2) are closely related to the functions, Qhe{R) 
and Fi^e{R) that we want to compute. When partial integration is applied to the variable 
Zi, boundary terms will appear since the integration area A covers only a part of the full 
droplet. The boundary terms can be expressed in terms of Q{R) and F{Ri, R2) and their 
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first two derivatives and give the relations 



Q^fl) = QiR) + c,{K)^^ + c,{R)^^^, (36) 
FUR) = F(fl,iS)+c.(fl)g+^^^^^ 




dRj dRl) ^ ' \dR1R2y ^ 



I d^F d^F \ I &^F \ 

+ci(^)c2(i?)b^ + 7^^ +^2(^) 77^ • (37) 



dR\R2 dRxR\) ^ '\dR\Rl^ 



R 



The subscript R means that both Ri and R2 should be set equal to R after differentiating. 
Eqs. (|3^) and (|37D are exact equalities, and the derivations are shown in appendix 0. We 
have used the notation 

The numerical method for obtaining Qhe{R) and Fi^e{R) is as follows: First gener- 
ate coordinates according to the absolute value \pLe\, and estimate Q{R) and F(-Ri,i?2) 
simultaneously in the way described above, but now with the advantage of having a gener- 
ating function that is symmetric in the electron coordinates. The computation is done for 
R, Ri, R2 = hk where h is a fixed grid spacing and k = 0,1, ■ ■ ■ , K. The maximum value 
K is taken large enough that the radii take values larger than the radius of the electron 
droplet, which for N electrons at filling fraction 1/m is y/mN. From the resulting one- 
and two-dimensional grid of data we estimate the derivatives by the formula 

^ ^^^^^ + ~ ~ (^^^ 

with similar expressions for higher derivatives. The values for QLe(-R) and Fi^dR) are 
found by use of Eqs. (^) and (pTf), and the value of AQi^e{R) by use of Eq. (|Tl]). In the 
computation we have used the two step sizes h = 0.1 and h = 0.2, and as will be discussed 
in detail in section 4, we have sufficient control over the errors introduced by the discrete 
differentiation. 
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3.4 The Jain quasi-electron 



We have also computed the charge and charge fluctuations for the wave function deflned 
by Jain. The wave function has the form [El 



Je 



k<l 

N / 1 



1 


7* 
Z2 

1 




1 








Z2 




Zn 


n 

k<l 


Zk — Zi)"^ 


N-2 

1 


Z2 




Zn 






{Zk - 


- ziT- 


N 




y 


n (^^ - 



j=l \ k<l:k,l=/=j 
'V/ k<l 



(40) 

(41) 
(42) 



Here V means projection onto the lowest Landau level. The importance of this projection 
1l5| , and it turns out that even without V most of the wave function 

Our calculations, which are performed for the 



has been studied 



is already in the lowest Landau level, 
unprojected as well as the projected state, give results that are in accordance with this. 

Following the general scheme for projection onto the lowest Landau level |T^, we let 
z* with the differentiation operator acting on everything except the exponential 

factor e"2 1]i=i l^'l^. This yields 



Je 



Zl 



k<l 



N 

rE 



ill 



i+3 



Zj Zi 



-E 



1 



+ 



m 



{Zj — Zk) J2kytj{Zj — Zk) 



(43) 

This expression is symmetric in all coordinates. In addition, the squared absolute value is 
the true electron density in the system, so in order to flnd Qje{R) and AQj^i^R) we can 
simply adopt the method used for the Laughlin ground state and quasi- hole state, Eq. ([T7|), 
now with pje = |V'JeP/-^je ^s the probability distribution used in the Metropolis algorithm. 
For the unprojected state the method is of course similar, and the expression in Eq. (|4^ 
now without the projection operator V, is used for pje- 



4 Results of charge mean value and charge fluctuation 
computations 

This section presents the results following from the calculations described in the previous 
section. For m = 3, i.e. for flUing fraction 1/3 of the lowest Landau level, we have 
calculated charge and charge fluctuations in the cases of a Laughlin quasi-hole, a Laughlin 
quasi-electron and a Jain quasi-electron. In all three cases the quasi-particle has been 
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located at the center of a circular electron droplet. For the quasi- hole we have made 
calculations for electron droplets consisting of 50 and 100 particles. Systems with 50, 
100 and 200 electrons were considered in the case of a Laughlin quasi-electron, whereas 
we did calculations for 50 electrons with the Jain quasi-electron in the system. We find 
that for all quasi-particles the numerical calculations reproduce well the expected bulk 
values of the charge mean values. There are finite-size effects, due to the limited number 
of electrons, which are most dominant for 50 electrons in the case of the Laughlin quasi- 
electron. However, these effects are constrained to the regions close to the position of the 
quasi-particle and to the edge of the electron droplet. For 100 and 200 electrons (and even 
for 50 electrons in the case of a quasi-hole or a Jain quasi-particle) the bulk values, —1/3 
for quasi-holes and +1/3 for quasi-electrons, are reproduced within small statistical errors. 
Charge is then measured in units of the electron charge. For the quasi-hole and Jain quasi- 
electron the calculations reproduce, again within small statistical errors, the expected bulk 
value zero for the charge fluctuations. However, for the Laughlin quasi-electron we find 
larger fluctuations than in the other two cases. At an absolute scale they are small, but 
they are significantly different from zero within the small errors of the calculation. 




Figure 3: Quasi-hole charge Ch(-R) for systems with 50 and 100 electrons. The horizontal line is —1/3. 

Fig. ^ shows the -R-dependence of the quasi-hole charge Ch(-R), for system sizes of 50 
and 100 electrons. The figure shows that there are three distinct regions. For small R there 
is a region where the charge builds up when R increases, which we clearly may identify with 
the location of the quasi-particle. There is an intermediate region where the charge seems 
to stabilize at a constant value, which is consistent with the expected bulk value —1/3 for 
the quasi-hole charge. For larger R there is a region where the charge again decreases to 
zero, and this we identify with the edge region of the droplet. The charge profile of the 
quasi-hole, for small i?, is essentially identical for 50 and 100 electrons, and that is also 
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(a) 



(b) 



Figure 4: (a): Quasi-hole charge fluctuations -Dh(i?) for systems with 50 and 100 electrons, (b): Enlarged 
vertical axis emphasizing the bulk behaviour of Dh{R) for 100 electrons. 



the case for the charge profile at the edge. The main difference between the two cases is 
the size of the intermediate region, the bulk region of the electron droplet. 

The distinction between the three regions for m = 3 is similar to what is seen in 
Fig. I^a for m = 1 . The main difference is the presence of oscillations for m = 3, in the 
charge profile for small and large R, which to some extent extends into the intermediate 
region. These oscillations are of the same form as has previously been found in numerical 
calculations of the charge density of Laughlin quasi-particles |]T3| . 

In Fig. ^ the charge fluctuations -Dh(-R) of the Laughlin quasi- hole are shown. Figure 

(a) compares the results obtained for systems with 50 and 100 electrons, whereas figure 

(b) shows an enlarged picture of the 100 electron case. The irregularities of the curve seen 
here are presumably due to the statistical fluctuations in the Monte-Carlo calculations 
and they give an indication of the size of these statistical errors. The results shown in 
Fig. I confirms the picture that the effect of the quasi-hole is restricted to a limited region 
around the origin. In the bulk of the electron droplet the charge fluctuations vanish within 
small statistical errors, whereas there are substantial fluctuations in the charge inside the 
quasi-hole and at the edge. According to the discussion in section ^, this is consistent with 
the assumption that the charge of the Laughlin quasi-hole is a sharply defined quantum 
number. 

The results presented so far are obtained for a discrete set of i?-values, R = hk with 
h = 0.2 and k = 0,1, ■ ■ ■ , K. The maximum number K is chosen such that hK is larger than 
the radius y/rnN of the electron droplet. Recall that lengths are measured in units of a/2£_b, 
with is SiS the magnetic length. The ground state data are for the case of 50 electrons 
obtained from 45 million electron configurations, and for 100 electrons from 17 million 
configurations. All quasi-hole data are found from 10 million configurations. For each data 
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set we have numerically estimated the standard deviation of the considered mean value. 
This enables us to estimate the numerical errors in the differences Ch(-R) = Qh{R) —Qo{R) 
and -Dh(-R) = AQh(-R) — ^Qo{R) as well. Table |^ shows some of the expectation values 
along with the computed statistical errors in the case of 50 electrons. The results are 
consistent with a quasi- hole charge that has the bulk value —1/3 and vanishing charge 
fluctuations. The estimated values for the standard deviation are somewhat larger than 
the irregularities of the plotted curve shown in Fig. ^ 



R 


Ch(i?) 




4.0 


-0.3219 ±0.0003 


-0.0028 ±0.002 


6.0 


-0.3334 ±0.0004 


±0.0003 ±0.005 


8.0 


-0.3371 ±0.0005 


±0.0000 ± 0.008 


10.0 


-0.3345 ±0.0005 


-0.0044 ±0.013 



Table 1: Charge Ch and charge fluctuations Z3h for the Laughlin quasi- hole in the case of 50 electrons. 
Both mean values and estimated errors are listed. The quantities Qh and AQh were obtained from 10 
million electron configurations, while 45 million configurations were used to find Qq and AQo- 
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0.3 
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Figure 5: Charge Clc(^) of a Laughlin quasi-electron. For < i? < 2.4 the curve is obtained by using 
Method 1 of section 3.3 for a system with 50 electrons. For R > 2.4 Method 2 was used, and the curves 
are for 50, 100 and 200 electrons. The horizontal line is 1/3. 



We will now turn to the case of a Laughlin quasi-electron located at the origin, and we 
consider first Fig. |] where the charge CLe(-R) is displayed. For R < 2.4 the curve is found 
for a system with 50 electrons using Method 1 of the previous section. The three different 
curves for R > 2.4 are for 50, 100 and 200 electrons and are found using Method 2. We 
observe again that there is a well defined region where the value of the charge is almost 
constant and agrees with the expected bulk value of 1/3 for the integrated quasi-electron 
charge. A more detailed presentation of the results show that for 100 and 200 electrons 
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this is indeed the case within the small errors of the calculation (to be discussed below) 
whereas for 50 electrons there is a small deviation which we assume to be a finite-size effect 
(see Fig. |^). Comparing with Fig. ^ we see that the curves for CLe(-R) are, apart from a 
sign, similar to the curves obtained for the quasi-hole. 




-0.002 



(a) (b) 



Figure 6: Laughlin quasi-electron charge fluctuations Dlc(^)- The single curve for < i? < 2.4 is found 



by using Method 1 of section 3.3 in a system with 50 electrons. The three curves for R > 2.4 are for 50, 



100 and 200 electrons, and Method 2 was used. 

Fig. H presents the calculated values of the charge fluctuations DLe(-R) of the Laughlin 
quasi-electron. (Also here the single curve shown for i? < 2.4 is obtained for a system of 
50 electrons with the use of Method 1, and for R > 2.4 the three curves are obtained for 
systems of 50, 100 and 200 electrons with the use of Method 2.) Figure (b) is an enlarged 
version of figure (a), now without the small R dependence. We observe that for the system 
sizes considered here, there are surviving fluctuations for the whole range of i?-values inside 
the electron droplet. This result is seen most clearly in Fig. ^d, where the irregularities 
indicate the size of the statistical fluctuations in the computation. A further discussion of 
uncertainties in the result will follow below. 

The behaviour of the charge fluctuation curve for the Laughlin quasi-electron is clearly 
different from that of the quasi-hole, which is shown in Fig. ^. For a system of 100 
electrons the quasi-hole charge fluctuations were seen to vanish in the bulk, within the 
small statistical errors of the calculations. A comparison between Figs, and with 
compensation for the difference in vertical scale, emphasizes the difference between the 
two cases. We thus conclude that there are charge fluctuations for the Laughlin quasi- 
electron that are larger and extend much further out than they do for the quasi-hole. Even 
though the expectation value of the quasi-electron charge is well defined and agrees with 
the expected bulk value (within the statistical uncertainty) in an intermediate interval of 
R, there are small but non-vanishing charge fluctuations which persist also here. In this 
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(a) DLe{R). 



(b) CLe(i?). 



Figure 7: Four statistically independent curves for the charge fluctuations Di^dR), (a), and the charge 
CLe(-R), (b) of a Laughlin quasi-electron. Each curve is based on approximately 48 million configurations 
for the quasi-electron data. The same ground state data, found from 125 million configurations, were used 
in all four cases. The lattice constant is h = 0.1 

sense the charge is not a sharply defined quantum number for the quasi-electron. However, 
there is a slow decay in the fluctuation curve for increasing R, and we cannot rule out 
that for even larger systems, with electron numbers N > 200, the value of the charge 
fluctuation will settle at the value zero further out in the bulk of the electron droplet. 
If that is the case, the quasi-electron charge is a sharp quantum number, but the size of 
the quasi-electron, as measured through these fluctuations, will be much larger than the 
size estimated from the charge expectation value or found by comparison with the charge 
fluctuations of the quasi-hole. 

The curves shown here for i? > 2.4 are based on the following number of electron con- 
flgurations: To flnd the quantities Qhe{R) and AQi^^R) we used 191 million conflgurations 
in the case of 50 electrons, 42 million for 100 electrons and 102 million for 200 electrons. 
The ground state quantities Qo{R) and AQq{R) were obtained from 125 million conflgu- 
rations for 50 electrons, 17 million for 100 electrons and 86 million electron conflgurations 
in the case of 200 electrons. For R < 2.4 we used 4.4-10^ electron conflgurations to flnd 
the quasi-clcctron quantities Qhe{R) and AQi^e{R). The calculations were done for values 
of R = hk, with k = 0,1, ■ ■ ■ , K , and h = 0.1 in the 50 electron case and h = 0.2 for the 
two other system sizes. We would also like to mention that in the case of 50 electrons, for 
a small range of R around R — 2.4, both the calculation methods. Method 1 and Method 
2, converged well and gave coinciding results for CLe(-R) and Di,e{R)- 

The conclusion for the charge fluctTiations suggested above is based on the assumption 
that the irregularities seen in the fluctuation curves give a measure of the errors introduced 
in the computation. However, in order to establish the conclusion more flrmly, we have 
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Figure 8: Charge fluctuations DLe(-R), (a), and charge Clc(-R), (b), found from two different lattice 
spacings h = 0.1 and h = 0.2 in the case of 50 electrons. 

considered numerical errors more carefully. The errors include accidental errors due to 
statistical fluctuations in the numerically obtained mean values, as well as systematical 
errors due to the finite lattice spacing h used in the numerical derivatives. Since the charge 
fluctuations, as well as the charge expectation values, have been calculated by a different 
method, which involves "negative probabilities" , the standard deviations of the calculated 
mean values are not so easily established, as in the case of the quasi-hole. In the case of 
the Laughlin quasi-electron the statistical errors therefore have been estimated simply by 
comparing the results of independent runs of the Monte-Carlo routine. To investigate the 
possible systematical errors we compare the results obtained by the use of two different 
step lengths and we compare the results of numerical differentiation with exact results, by 
applying the methods to a representative test function which can be handled analytically. 

The estimate for the statistical errors is obtained by dividing the numerical data for 
Q{R) and F(i?i,i?2) into four statistically independent piles, each consisting of approx- 
imately 48 million electron configurations. For each pile we have computed CLe(-R) and 
Di^e{R), using the same ground state data in all cases. Fig. ^ shows the four fluctua- 
tion curves found in this way for the case of 50 electrons, while Fig. ^ shows the four 
independent charge curves. For a given R we estimate the errors in the mean values as 
half the difference (-DLe(-R))max — (-DLe(-R))min5 and similarly for Ci^^{R). For the charge 
fluctuations this number is from the figure seen to typically take the value 0.0005, and 
implies for example for R = 7 that the mean value found is approximately 10 times larger 
than this error. Fig. ^ also shows that the estimation of the errors obtained through the 
independent runs gives essentially the same results as the estimation obtained visually from 
the irregularities of a single curve. As judged from Fig. ^ the statistical errors for larger R, 
in the 100 and 200 electron curves, may be slightly larger than displayed by the 50 electron 
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curves. However, the value of the fluctuations still are significantly different from zero in 
most of the electron droplet. 

A similarly estimated error in the charge expectation value CLe(-R) is approximately 
0.001. A detailed presentation (not included here) of the results in the cases of 100 and 200 
electrons shows that the calculated charge expectation values do not deviate significantly 
(within this error) from 1/3 in the bulk of the electron droplet. 




{8i)F{Ri,R2) (b) /i(i?i,i?2) 



Figure 9: Comparison between (a): The Monte Carlo estimated function F{Ri,R2) for 50 electrons 
and lattice spacing h — 0.2, and (b): The analytically found test function /i(i?i,i?2) from Eq. (|4^), also 
calculated for 50 electrons. In both cases the factor N{N — 1) has been divided out. 

We will now discuss the systematical errors introduced by the finite lattice spacing h. 
For small h the lattice spacing will give errors in the numerical derivatives, Eqs. ( ^6|) and 
(^), which are proportional to h^. The above results for 50 electrons are all found using 
h = 0.1. In addition we have performed computations with h = 0.2, and Fig. | compares 
curves for Di^^iR) and CLe(-R), obtained with these two values of h. We observe that the 
differences are similar in size to the statistical error that we have already considered. Due 
to the /i-dependence of the error for small h we expect that a further reduction in the 
lattice spacing will introduce differences which are smaller than the differences between 
the results obtained with h = 0.2 and h = 0.1. Results for h = 0.02 seem to confirm this. 

To make an independent check of the error introduced by the discrete differentiation, 
we have applied this operation to a test function with a similar form as F{Ri,R2). From 
Eqs. (|32D, (p^) and (0) we observe that, apart from the factor nili(kiP — 1) the function 
F(-Ri,i?2) is nothing else than the integrated two-particle density of the ground state for 
filling fraction 1/3. We therefore expect that except for a scaling in distance, the function 

m„R2) = ^^^~^^ I d'zj dhJd'<^''-'h\H', (44) 

il JAi JAi J 



20 



will mimic the main properties of the function F(i?i,i?2) reasonably well. Here ipi is the 
wave function for the m = 1 ground state (0), and Ji is the corresponding normalization 
integral. That this anticipation is indeed fulfilled is seen in Fig. where the analytically 
determined function /i(i?i,i?2) is compared to the numerically found function F{Ri,R2). 
Both cases are for 50 electrons, and in each figure the surface starts from on the Ri 
and R2 axes, and then smoothly builds up to the value N{N — 1) as the edge of the 
electron droplet is reached. For the test function /i(i?i,i?2) we have then computed the 
right hand side of Eq. ( pT]) in two ways: We have computed it analytically, and in addition 
we have calculated it by use of the discrete differentiation, Eq. (^9]), with lattice spacing 
h = 0.1. We found errors in the result obtained with the discrete derivatives which were 
approximately 0.0005 in the central part of the droplet. This agrees well with the results 
found from runs of the Monte Carlo routine with the two different spacings, h = 0.2 and 
h = 0.1, and are similar to the estimated value for the statistical errors. 

Based on these estimates of the errors introduced in the computation we find it rea- 
sonable to conclude that the deviations from zero seen in Fig. ^ for the charge fluctuations 
are real and not due to the errors introduced by the calculations. 




(a) (b) 

Figure 10: Charge Cjo(i?) of a Jain quasi-electron in a system with 50 electrons, (a): Projected state, 
Eq. (H). (b): Unprojected state, Eq. (|2|). 

As noted above, we have also studied the charge and charge fluctuations for the quasi- 
electron defined by Jain's wave function. This investigation was prompted by the behaviour 
found for the Laughlin quasi-electron, as a wish to compare the different proposals for the 
wave function. For a Jain quasi-electron located at the origin calculations have been 
performed for the projected state, Eq. (|43|), as well as for the unprojected state, Eq. (|4^). 
In both cases we have considered only a system of 50 electrons. 
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(a) 



(b) 



Figure 11: Charge fluctuations Djo(i?) of a Jain quasi-electron in a system with 50 electrons, (a): 
Projected state, (b): Unprojected state. 



Fig. compares the charge Cje(-R), as found for the projected state, to the corre- 
sponding quantity for the unprojected state, with the results for the projected one in the 
figure (a). We readily notice that there is almost no difference between the two figures, 
i.e. projection onto the lowest Landau level does not make any important difference in the 
present context. In both cases there is a well defined intermediate region where the curve 
is consistent with the expected integrated charge value 1/3 (times the electron charge). 
When compared to Fig. |^ we notice that the amplitude of the peak for small R is larger 
for the Jain quasi-electrons than it is for the Laughlin quasi-electrons, but apart from this 
the differences in charge expectation values are rather small for the two definitions of the 
quasi-electron wave function. 

This similarity changes when we consider the charge fluctuations. For the Jain quasi- 
electron these are displayed in Fig. |lT], and Fig. shows the same results with the vertical 
axes enlarged. We observe again that the projection onto the lowest Landau level has a very 
small effect on the result. But the most compelling observation is that there exists a small 
region, which we interpret as corresponding to the bulk of the droplet, where the charge 
fluctuations are very small, indeed consistent with the value zero within small statistical 
errors (see the discussion below). The range of R for which this is true is about the same 
as the range where the charge mean value, shown in Fig. ITDl is close to the bulk value 1/3. 
This behaviour is similar to the results we found for the Laughlin quasi-hole, and differs 
from what we saw for the Laughlin quasi-electron. 

The results presented here are based on the following set of data: For the ground state 
we used 45 million electron configurations, whereas we for the Jain quasi-electrons used 77 
million configurations for the projected state and 34 million for the unprojected. We have 
performed calculations for R = hk with h = 0.2 and k = 0,1 ■ ■ ■ , K. For each set of data 
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(a) (b) 



Figure 12: Fluctuations -Djo(-R) of a Jain quasi-electron in a system with 50 electrons, (a): Projected 
state, (b): Unprojected state. 



R 


CUR) 




4.0 


0.3640 ±0.0002 


-0.001 ±0.001 


6.0 


0.3281 ±0.0002 


±0.001 ±0.003 


8.0 


0.3355 ±0.0003 


-0.001 ±0.005 


10.0 


0.3383 ±0.0003 


±0.005 ±0.007 



Table 2: Charge Cjo and charge fluctuations Djc for the Jain quasi-electron in the case of 50 electrons. 
Both mean values and estimated errors are listed. 

we estimated the numerical standard deviations in Cje and Dj^,. Some of the results are 
listed in Table 0. The table shows, for the listed values of R, that the calculated values for 
both Cjo and Dje are consistent with the expected bulk values within the estimated errors. 

5 Statistics parameter of the Jain quasi-electrons 

The studies of the charge fluctuations presented in this paper have been motivated by the 
results of Ref. where the charge and statistics parameters of the Laughlin quasi-electrons 
were studied numerically, by Berry phase calculations, for systems of up to 200 electrons. 
Whereas the expected charge value, in those calculations, was well reproduced, no stable 
value for the statistics parameter was found. This was clearly different from calculations of 
the statistics parameter of the quasi-hole. Our hypothesis has been that the discrepancy 
between the quasi-hole and quasi-electron results is due to more long range fluctuations 
in the correlation functions for the latter. This conclusion seems to be supported by the 
results of the charge fluctuation calculations presented in this paper. 



23 



The alternative definition of tlie quasi-electron wave function, given by Jain, seems to 
be better behaved as far as the charge fluctuations are concerned. Although we have not 
examined this wave function in an equally detailed way as the one introduced by Laughlin, 
this seems to be a reasonable conclusion, based on the results of the calculations on the 
50 electron system. Thus, the fluctuations relative to those of the ground state, when 
moving away from the center of the quasi-electron, seem to be more rapidly damped for 
Jain's wave function than for the wave function defined by Laughlin. To pursue this point 
further we will now present the results of numerical calculations of the charge and statistics 
parameter, as extracted from Berry phases, in the case of Jain's quasi-electron. The results 
presented are for a system of 74 electrons at filling fraction 1/3, and for simplicity we have 
used the unprojected quasi-electron wave function. We will show that in this case a rather 
well-defined value is found both for the charge and statistics parameters. The value of the 
charge parameter is close to the expected value 1/3, and the statistics parameter is close 
to —1/3. The sign of the latter is the opposite of that of the quasi-hole and it is also the 
opposite of what is expected from a simple model of the quasi-electron as a charge-fiux 
composite. 

Let us briefiy review how charge and statistics parameters are extracted from Berry 
phases. The idea, as originally set forth by Arovas et.al. in Ref. 0, is to let the position 
parameter zo that labels the quasi-particle state traverse a loop in the plane. The Berry 
phase 0] associated with this motion can then be computed. This phase is in turn 



interpreted as an Aharonov-Bohm phase for the unknown charge q of the quasi-particle 



encircling the known magnetic fiux, and the charge is extracted. Thus, there is a non-trivial 
interpretation that lies under this way of determining the charge. To find the statistics 
parameter one computes the Berry phase associated with two quasi-particles encircling 
one another, and interprets the two-particle contribution to the Berry phase as an anyon 
interchange parameter. 

To be more specific, let \zq) be the normalized state corresponding to a single quasi- 
particle located at the position zq, and suppose the particle is moved around a closed loop 
parameterized by zq = re^'^ with running from to 271. The Berry connection is then 
defined by 

A,{r)=t{zo\d^\zo), (45) 

and the Berry phase is the integral of the Berry connection along the path. This phase we 
relate to the Aharonov-Bohm phase for a charged particle in a (uniform) magnetic field. 
If the path is left handed relative to the direction of the magnetic field, as is the case in 
our calculations, the Aharonov-Bohm phase associated with the propagation of the charge 
around the loop is 

7 = 27iqr^. (46) 

Charge is here measured in units of the electron charge, and r is the dimensionless radius, 
measured in units of v^^b, of the circular loop. The charge of the quasi-particle is now 
determined by setting the Berry phase equal to the Aharonov-Bohm phase. If the Berry 
connection depends on r = \zo\ but not on (due to rotational invariance), the Berry 
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phase is (3i{r) = 27iAi{r) and the charge parameter is 

,(r) = i|<4. (47) 

In general this definition will give an r-dependent charge for the quasi-particle, but far 
from the edge of the electron droplet g(r) is expected to settle at a constant value. 

Suppose there are two quasi-particles in the system, with positions ±2:0, so that the 
parameterization zq = re**^, with now running from to vr, describes a counterclockwise 
interchange of the two particles. We can define the Berry connection associated with this 
interchange as 

A2 = i{zo,-zo\d^\zo,-zo), (48) 

where \zo, —zq) is the two quasi-particle state. Subtracting the single-particle contributions 
we have that u, defined as 

u{r) = -{P2{r)-2p,{r)), (49) 

TT 

with as the integral of A2, can be identified (for large separation r) with the anyon 



statistics parameter of the particles |23- Eqs. ( ^7|) and (|^) constitute the basis for the 
discussion below. 

Before proceeding, we would like to make a comment on the sign convention used in 
the definition of the statistics parameter (|^, since the sign is of some importance for the 
discussion. The sign is fixed by orienting the loop used for the Berry phase calculations 
in a positive direction relative to the product eB of the electron charge and the external 
magnetic field. Thus, it is independent of the sign of the charge qe of the quasi-particle. 
However, another convention is possible, and maybe even more useful. If the orientation of 
the loop in the calculation of the statistics parameter is fixed relative to qeB rather than 
eB, then the sign of the statistics parameter has a direct physical significance. We may 
write the parameter with this new convention as 

The sign of ui is then determined by the relative sign between one-particle and the two- 
particle contributions to ^2- If ^1 is positive there is an effective repulsive (statistical) 
interaction between the two quasi-particles, and if it is negative there is an attractive 
interaction^. 

^For anyons in the lowest Landau level the parameter is normally restricted to the interval < Vi < 2. 
However, there is a natural extension of this to the interval — 1 < j^i < cxd. Negative values then corresponds 
to anyons with an additional attraction which gives a singular (but normalizable) short range behaviour 
of the wave function. If vi is larger than 1 there is a repulsion corresponding to the exclusion of one or 
more of the lowest relative angular momentum states. (The Laughlin states of filling fraction 1/m are of 
this kind.) The statistics parameter vi with the given sign convention is identical (up to a constant shift) 
to the one-dimensional statistics parameter introduced in terms of algebraic relations between observables 
of the system . It is also identical to the exclusion statistics parameter defined by state counting in the 
many quasi-particle space |2l[| . 
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A Jain quasi-electron located at the origin is described by the wave function of Eq. (|^). 
To find the charge and statistics parameters we need to translate the quasi-electron to the 
position zq without moving the circular electron droplet itself. Since the expression for such 
a translated quasi-electron (to the best of our knowledge) does not exist in the literature, 
and since it does not follow trivially from (^) , we will include a discussion of how to solve 
this problem. 

Notice that apart from the projection operator V, the wave function (|40|) for m = 1 
describes a filled lowest Landau level with a single electron pushed up to the next Landau 
level, where it occupies the single-electron coherent state localized at the origin . If this 
single-electron state is translated to the point zo it is described by 



fzoiz,Z*) 



1 



■^{zz*+zoz*) 



TC 



(51) 



Here z is the electron coordinate, and the wave function is normalized to unity in the state 
space of a single electron. The coherent states are known to be maximally localized, so if 
one of the considered electrons now occupy fzo{z, z*) the m = 1 Jain wave function will 
have an excess charge equal to the electron charge accumulated close to the position zq. 
Hence the quasi-electron has been moved from the origin to Zq. 

For m 1 the situation is a bit more complicated. The VanderMonde determinant 



A = Ylizk-zi] 

k<l 



(52) 



raised to the power (m — 1), which is the new factor in the wave function as compared to 
the m = 1 function, will push the electrons apart. The detailed effect of this pushing is not 
completely known, but we notice that A treats all electrons symmetrically. This implies 
that a quasi-electron that for m = 1 was located at Zq now is moved, assuming that we 
keep the form in Eq. (^) of the coherent state. We do not have a simple argument to 
determine how the quasi-electron position will depend on m, but numerical studies show 
that the new position is close to mzo. This scaling with m of the quasi-electron position 
can be avoided if we replace zq with zo/m in Eq. (^Tj). Hence for a Jain quasi-electron 
localized at the position zq we use the wave function 
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(53) 



That this wave function indeed has excess charge localized around the position zq is shown 
in Fig. O. The figure shows a cut from the origin of the electron droplet and through the 
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point = 3, of the (non- normalized) single electron density^ 

with C as a constant. We see that the density profile, with its dip close to Zq and two 
peaks on each side of the dip mimics that of a coherent state in the first Landau level, 
Eq. (|51|), localized at zq = 3. This figure supports the definition we have used for the wave 
function describing a translated quasi-electron. 



Figure 13: Cut in the radial direction of the single electron density p{zi, z*), Eq. (|5^), for a Jain quasi- 
electron. The cut starts at the origin and goes through the parameter point zq = 3. System size is 50 
electrons. For computational simplicity the density is not normalized. 



For two quasi-electrons at the positions ±zo we use the wave function 



m— 1 



m ' 
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This function treats the two quasi-electrons symmetrically. 

Before we proceed to the numerical results we need to establish the relations between 
the desired Berry connections, given by Eqs. (^) and (^Uf), and the normalization integrals 
of the two wave functions above. The state \zq) in Eq. (^) generally can be expressed as 



K 



Vi fc=o 



(56) 



^The numerical method used to find p{zi, z^) will be describe below. 
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where \k) are orthonormal basis states, are expansion coefficients and 



K 

I = T.^^''W?: r=\zol (57) 

A;=0 

is the normahzation factor. The Berry connection can be related to / by the expression 

A(r) = r^-^ln/. (58) 
dr"' 

This relation relies on the property that I depends only on the absolute value r of and 
not its phase. Since an expansion of the wave function as a power series in Zq can be 
shown to be an expansion in terms of orthogonal total angular momentum eigenstates, in 
exactly the same way as in Eq. (|56|), the relation ([58|) holds, with 



Il = /d2^z|^L(^i,---,z^)p. (59) 



For the two quasi-electron state an expansion in Zq is again an expansion in terms of 
orthogonal angular momentum eigenstates. However, the expansion now has the form 

w-z,) = ^j:{z^r^w\k). (60) 

V J fc=o 

The lowest power of Zq is not in this expansion, which means that the wave function 
^/'jg contains an unphysical singularity for zq = 0. This singularity should be removed, and 
we can achieve this by using a complex normalization factor. This changes the relation 
between the Berry connection and the normalization integral. We find 

A,{r) = r'^^\nll-1, (61) 

where 

Il = ld''^z\^liz,,...,z*^)\' (62) 

is the usual normalization integral. To summarize, the charge and statistics parameter are 
related to the normalization integrals in Eqs. ( p9D and (p^ ) by 

^ = (63) 



d 



r 



2 ^ /i„ 7-2 



dr 



(ln/J,-21nJ]j-l. (64) 



We are now ready to present the numerical method used to find the normalization 
integrals Jjg and Jjg as a function of r = |zo|- Monte Carlo integration is used, and for the 
generating probability density p{zi, ■ ■ ■ , z]^) used in the Metropolis algorithm we have the 
following requirements: p should represent the integrands properly so that the numerical 
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uncertainty becomes as small as possible, and it should not depend on zq. To satisfy these 
criteria, we have used the squared ground state wave function as the probability density, 
that is 

p(^i,---,^;v) = ^e-Stil^»l>r- (65) 

Here Iq = J d^^2e~^i=i '^''^|Ap™. The wave function for a single quasi-electron, i.e. the 
function that enters the expression for Jj^, can be written as 

= e-^Sr^.l-I^A-E(4 - ^-H-iy'-'^ (66) 

i=l ^ ^ 

where Mj is the determinant arising when we remove the first row and the i'th column 
of the original determinant in Eq. ([53|). The wave function for two quasi-electrons can 
be similarly expanded, although now a double sum will appear for the sum in Eq. ([66|). 
The great numerical advantage we have achieved by our choice of probability density p is 
that for each specific electron configuration, the number of multiplications needed to find 
the ratio ^ is proportional to rather than A^^, the latter being typical for Mi and A 
separately. Here is the number of electrons. The reduction by the factor N dramatically 
reduces the required computer time. We generate electron configurations according to the 
probability density p above. For zq = r = hk, with h fixed and k = 0,1, ■ ■ ■ , K, we estimate 
/jg(r) and Ije{r). Notice that the two normalization integrals are found from the same set 
of electron configurations, and hence have correlated uncertainties. As noticed in Ref. 
this reduces the numerical error when we find u from the difference in Eq. (Q) since the 
numerical errors have a tendency to cancel. 

The numerical method used to compute the single electron density in Eq. ( |MD is basi- 
cally the same. However, in this case the coordinate zi as well as zq needs to be treated as a 
parameter in the numerical calculations. In this case we have therefore generated electron 
configurations by using the square of the (A^ — l)-electron ground state wave function as 
the probability density. 

The results of our Berry phase calculations will now be presented. For 74 electrons we 
have determined the charge q and statistics parameter u according to Eqs. ( |53D and (p^). 
The computations are done for m = 3, i.e. for 1/3 of a filled Landau level, and the required 
integrals l}e{f) and /je('") are computed for the parameter r taking the values r = hk, with 
h = 0.2 and k = 0,1, ■ ■ ■ , K. Fig. O displays the results. 



Fig. 0a shows the charge q as given by Eq. (pSj) as a function of r, the dimensionless 
distance from the origin to the quasi-electron. We see that q is (almost) constant and 
stays close to the plateau value 1/3 all the way from the origin and until close to the edge 
of the electron droplet. For 74 electrons the latter has a radius of 14.9. This constant 
behaviour means that the bulk value of the charge, as extracted from Berry phases, is 
well defined. However, the small deviation from the expected value 1/3 is found to be 
statistically significant. It is interesting to notice that a similar deviation was seen for the 
charge of the Laughlin quasi-electron in Ref. . In that case it was interpreted as a finite 
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size effect that would vanish when N oo, since the deviation was seen to become smaller 
as was raised. For the present case, a calculation for 50 electrons gave the same value 
of the plateau as does Fig. |H|a, hence we have no indication that the deviation is due to 
the finite size of the electron droplet. We are aware that it may be related to the specific 
way we have defined the wave function for the translated quasi-electron in Eq. (|53|) , but 
we have not studied this point thoroughly. 




Figure 14: (a): The charge q{r) as a function of the quasi-electron's distance from the origin, (b): The 
statistics parameter iy(r) as a function of half the distance between the quasi-electrons. The system has 
74 electrons, and for comparison we have included the horizontal lines 1/3 and —1/3. 

Notice that even for r —>■ 0, q stays close to the plateau rather than dropping to zero 
as does the integrated single-electron density; compare to Fig. [1^. This emphasizes that 
the present definition of the charge does not necessarily give the same charge as the one 
found by integrating the single-electron density. When we interpret the Berry phase as 
an Aharonov-Bohm phase, the behaviour as r — > can in the present case be understood 
in the following way: The Aharonov-Bohm phase 7 associated with the loop zo = re*"^ is 
proportional to the encircled magnetic flux. The curve in Fig. 14a shows 2^7? which then 



is a constant. This holds even if the charge is not truly point like, because all parts of the 
smeared out charge is moved around in equal loops. 



We now turn to Fig. |T4|b where the statistics parameter i^(r), as given by Eq. (|6^), 
is shown as a function of r. The figure shows that there is a range of r for which the 
value of u is approximately constant. That is, we see rather clear signs of a well defined 
bulk value for the statistics parameter of the Jain quasi-electrons. To give unquestionable 
evidence for the existence of such a plateau, and also to determine the precise value of the 
parameter, would require more detailed studies than we have done here. However, we do 
observe that the value —1/3 seems likely to be the value of the statistics parameter u. This 
equals minus the statistics parameter of the Laughlin quasi-holes. Let us mention that we 
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(a) Statistics parameter for Jain quasi- 
electrons. 



(b) Statistics parameter for Laughlin quasi- 
electrons. 



Figure 15: Comparison between curves that define the statistics parameters of (a): Jain quasi-electrons 
(|64|), and (b): Laughhn quasi-electrons. For the Jain quasi-electrons we have used a system with 74 
electrons, whereas the Laughlin quasi-electron curves are for systems with 20, 50, 75, 100 and 200 electrons. 
Of these the 75 electron curve is farthest to the right, and the figure is from Ref. 



also for a system with 50 electrons have seen clear indications of this. 

The results seen here for the Jain quasi-electrons are similar to the results obtained in 
Ref. IP for the Laughlin quasi-hole. But the most striking observation is the big difference 
between the results found in this reference for the Laughlin quasi-electron, and the present 



results for the Jain quasi-electron. Fig. 15 shows a comparison between the curves that 



are used to define the statistics parameter, as found from Berry phases, for these two 
definitions of the quasi-electron wave function, with Jain's definition used in figure (a) and 
Laughlin's definition used in (b). The latter figure is from Ref. and shows results for 
the system sizes 20, 50, 75, 100 and 200 electrons. The 75 electron curve, which is most 
appropriate for comparison with our Jain quasi-electron curve for 74 electrons, is the one 



farthest to the right. The difference between Figs.|T5|a and p!5| b is striking; whereas the 
curves for Laughlin quasi-electron do not show signs of a well-defined statistics parameter, 
that is clearly the case for the Jain quasi-electron. 

The signs of the statistics parameters of the quasi-electrons deserve special attention. 
We note that they are different for Jain's and Laughlin's wave functions. Let us re-express 
the results in terms of the parameter ui which measures the repulsion or attraction between 
the quasi-particles. For the quasi-hole the value is z/i = +1/3, which corresponds to an 
effective repulsion between two quasi-holes. Compared to the quasi-hole the Jain quasi- 
electron has the opposite sign for the charge as well as for the statistics parameter u [Q. 
This means that ui is positive also for the Jain quasi-electrons. The opposite sign for the 
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statistics parameter, as may be indicated by the Berry phase calculations for the Laughlin 
quasi-electron, would correspond to an effective attraction between the quasi-electrons. In 
the anyon representation this attraction would be represented by a singular anyon wave 
function 0. 

A simple picture of the quasi-particles represented as charge-flux composites indicates 
that (the non-integer part of) the statistics parameter ui should have opposite sign for 
the quasi-hole and the quasi-electron. This is because the sign of both the flux and the 
charge is reversed for the quasi-electron compared to the quasi-hole. This gives the same 
value for the parameter u, but the opposite sign for ui, due to the different signs of the 
charges. This simple argument is supported by numerical studies of the quantum Hall 



effect on a sphere |^ . State counting gives an exclusion statistics parameter for the quasi- 
electrons which corresponds to z/i = 2 — 1/3, as compared to ui = +1/3 for the quasi-holes. 
The fractional part is negative, but it is compensated by a positive integer part to give 
altogether a repulsive parameter. It is of interest to note that the results obtained for Berry 
phase calculations of the statistics parameter do not fit this expected value of the physical 
quasi-electron, neither with Jain's nor Laughlin's wave functions for the quasi-electron. 



6 Conclusions 

To summarize, we have examined expectation values and fluctuations for the charge of 
quasi-holes and quasi-electrons in the quantum Hall system with filling fraction 1/3. The 
study has been motivated by the asymmetry between quasi-holes and quasi-electrons that 
has previously been found in Berry phase calculations of the statistics parameter 0. We 
have in particular been interested in examining to what extent the quasi-electron charge 
can be regarded as a sharply defined quantum number. To study these problems we have 
calculated charge fluctuations, as well as charge expectation values, measured relative to 
the ground state. We have done computations for a Laughlin quasi-hole, and for quasi- 
electrons both with Laughlin's and Jain's definition of the wave function. In all cases the 
quasi-particle has been located at the center of a circular electron droplet. 

For the quasi-hole charge we have found numerical results that are consistent with the 
expected bulk value —1/3 (times the electron charge) and charge fluctuations consistent 
with vanishing fluctuations in the bulk. There are effects due to the finite size of the quasi- 
hole and to the finite size of the droplet, but in an intermediate interval the deviations 
from the expected values vanish within small statistical errors. The computations have been 
done for systems with 50 and 100 electrons, and the results clearly confirm the conclusion 
that the quasi-hole charge is a sharp quantum number. 

For the Laughlin quasi-electron the calculations similarly reproduce the expected bulk 
value of the charge, which is 1/3 times the electron charge. But for the charge fluctuations 
of this excitation we find fluctuations which survive throughout the whole electron droplet. 
At an absolute scale the values are small, but we have shown that they are significantly 
different from zero. We have studied systems with 50, 100 and 200 electrons, and conclude 
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that results obtained for these system sizes do not unambiguously confirm the charge of 
the Laughlin quasi-electron to be a sharp quantum number. The deviation from zero which 
are found for the charge fluctuations may be a finite-size effect, but the fluctuations then 
are sufficiently long range that even a system of 200 electrons is not sufficiently large to 
see this clearly. 

For a system with 50 electrons we have also studied the charge and charge fluctuations 
of a Jain quasi-electron. In this case we again reproduce expected results in the bulk, 1/3 
for the charge and vanishing fluctuations. There are small deviations, but these are not 
significant, within small statistical errors. We also find that in this context there is almost 
no difference between the results obtained for the Jain wave function projected onto the 
lowest Landau level and the unprojected state. This confirms results presented elsewhere 
in the literature, which show that the unprojected state is almost entirely in the lowest 
Landau level. 

As a further examination of the relation between sharpness of the charge and a well 
defined statistics parameter, we have computed Berry phases for the Jain quasi-electron. 
This gives a supplement to the previously done calculations for the Laughlin quasi-electron 
0. We have considered unprojected states for a system with 50 electrons, and have ex- 
tracted charge and statistics parameters in the usual way. The charge parameter calculated 
in this way has a well defined bulk value with only a small deviation from the expected 
value 1/3. The small deviation may be due the definition used for the wave function of a 
quasi-electron translated to arbitrary position, but we have no firm conclusion about this. 

As the most interesting result of the Berry calculations, we found that the statistics 
parameter of the Jain quasi-electron is much more well behaved than the analogous quantity 
for the Laughlin quasi-electron. With the convention used here, our results clearly indicate 
that the value of this parameter (for large separations of the two quasi-electrons) is z/ = 
— 1/3, which is minus the statistics parameter of the quasi-holes. However, the sign of the 
statistics parameter does not agree with the expected one, and is not consistent with the 
value —2 + 1/3 indicated by state counting of numerically determined energy levels for 
quantum Hall states on a sphere. Thus, for neither of the suggested wave functions, due to 
Laughlin and Jain, the expected statistics parameter of the physical quasi-electrons seem 
to be correctly reproduced. 
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A Numerical techniques 



This appendix reviews important aspects of the Metropohs algorithm and Monte Carlo 
estimate. 

In section ^ we give several examples of functions used to generate electron configu- 
rations according to the Metropolis algorithm, e.g. Eq. (|l5l). We will briefly review the 
technical details of this method. Suppose then that the desired (real) probability distri- 
bution is given by p{zi, ■ ■ ■ , z^), and there is a given configuration {zi^a}iLi- To find the 
next configuration we loop through the electrons, and for each i we randomly choose a test 
coordinate = Zi^a + such that Azi and —Azi are equally probable. We then compute 
the ratio 

Ji = —( z ^ (67) 

If the number /j is larger than a randomly generated number between and 1, then we 
accept the test coordinate and set Zi^a+i = z\. Otherwise Zi^a+i = Zi^a- The procedure 
ensures that the configurations are generated according to the desired probability distribu- 
tion p because the principle of detailed balance is satisfied by the transition probabilities; 
the ratio of jumping from Za to Zi, or from z^, to Za equals the ratio p{zb) / p{za) ■ 

For the cases we have considered it is of special interest to notice that an overall normal- 
ization factor in p is irrelevant for the Metropolis algorithm since the latter only considers 
ratios between different value of p. In all our calculations of charge and charge fluctuations 
we have benefited from this in the sense that we have used probability distributions for 
which we did not know how to analytically find the normalization factor. For instance the 
factor Iq in Eq. (|TB[) is not known exactly. It was essential that we could do this, because 
the computer time needed to obtain well converged results highly depends on the choice 
of probability distribution. The time is reduced when p has a behaviour similar to the 
actual integrand of the problem under consideration, because the standard deviation of 



the Monte Carlo estimate is given by ^^JV{g/p)/n, with g as the specific integrand, V{g/p) 
as the variance of the function g/p, and n the number of Monte Carlo steps. 

B Relations between R2) and QLe{R), Fi^e{R) 



In this appendix we will show how to derive the expressions in Eqs. ( |36| , pT] ). 

We start by considering the relation between Q{R) and QheiR)- According to Eq. (^ 
the desired charge expectation value can be written 

QUR) = ^ f d\ fd'^''-'h\ijU^„---,z%)\' (68 



he J A 
rR 

2tt / drrh{r), (69) 
Jo 
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where we have defined 

h(r) = 

2tt Jo 



Hr) = ^ r dippizi,ziy, zi = re'^ (70) 



with 



p{zi,zl) = N^e-\'^\'d,At (71) 



'Le 

N 



J d'z, ■ ■■d'zr.e-^l^^'^^'' l[{\z,\' - 1) n - zj^. (72) 



k=2 i<j 



It is important to notice the ranges of the counting variables in this expression, which is 
obtained using integration by parts for all coordinates that are integrated over the entire 
complex plane. The function p{zi,zl) is the exact single-electron density. We now define 
the quantities 



1 r N ^ 

piz^^zl) = N— / d'z, ■ ..d'z^e-T...^\^'\' n(k.P - 1) n 1^^ - ^l'"'' (73) 

-'Le •' k=l i<j 



2tx Jo 

and 



h{r) = ^ r difpiz^, zl) (74) 



QLeiR) = 27r / drrh(r). (75) 
Jo 

The latter quantity is identical to the one defined in Eq. (pi]). Comparing Eqs. (^) and 
([f3|) we observe that the relation between p and p can be written 

p{z„zl) = e-\^^\'d^Al (el-l^^p-i-^p(zi,^t)) • (76) 
It is then straight forward to show that the radial functions h and h are related by 

Hr) = Mr) + 1 A ( ^h(r) + . (77) 

r dr — 1 4 dr — 1/ / 

This implies that 

QUR) = 27r ["^drrhir) (78) 
Jo 



2n l''drrhir) + 2TT(-^hiR) + \R(^j^^^] V (79) 



r=R/ 
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since the contribution from r = is zero. Notice that the appearant singularity at r = 1 is 
artificial since the factor — 1 in the denominators is canceled by the same factor in the 
numerators; recall the definition of h{r). Using the fact that 

we finally obtain the advertised relation 

QUR) = QiR) + ^i(R)^^ + (81) 

with 

To find the relation between Fi^e{R) and F(i?i,i?2) we perform a calculation that is 
similar in spirit, although technically more complicated than the one above. First we define 
the radial function 

^2(ri,r2) = r dif,dip2P2{zi, Z2, zl z;); Zi = re''^\ Z2 = re^^^ (83) 



47r2 JO 

with 



P2[z,,Z2,zlz;) = N{N-l)^e-\''\"-\''\"d,Ald.A*, (84) 
X / d'z,---d'z^e-^^-s\^^\' l[i\z,\' -l)l[\z,- z/"' (85) 

A:=3 i<j 



as the true two-particle distribution function. We can then write the desired function in 
Eq. (H) as 

FUR)=4:7c^f dn f dr2nr2h2{n,r2). (86) 

JO JO 

Let us now define the auxiliary quantities 

P2{z^,Z2,zl,z;) = iV(iV-l)-L (87) 

N 

X I d^z^---d^ZNe ^«=i'-' Y\(\zkf -l)Y\\zi- Zjl"^"", 



J d'zs ■ --d'z^e-^i^ i-i^ n(k^r - 1) n - 

k=l i<j 

/i2(n,r2) = dipid(p2P2{zi,Z2,zl,z^), (89) 

and 

F{Ri,R2) =47v^ dnn / dr2r2/i2(ri,r2). (90) 

JO JO 
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Analogous to what we saw above, also now the relation between the two-electron density 
P2 and the quantity p2 is easily established, and it reads 

A similar exact relation exists between the radial quantities h2{ri,r2) and h2{ri,r2). We 
rewrite it into a form suited for integration and find 

h2{ri,r2) = h2{ri,r2) 

,1/2 ^ dg2\ 1,2 "9 / dg2\ 



4 ri dri \ dri J 4: r2 dr2 \ dr 



2 ) 

11 d ( d d , 2 \ 11 d ( d d , 2 \ 

+ 7 ^^'^'^'^'^292) + 7 r-2 {r^g2) 

4 rir2 OTi \ OTi or2 J 4 rir2 o'r2 \ or\ or2 J 

1 "9 <9 / 2 2 \ 
rir2 OTi ^2 

1 1 a as d 

+ ^^l o o ^2T^^2- (92) 

lo rir2 ari ori or2 or2 

We have in this expression used the notation 5'2('"i5 ''^2) = / a_i w :a_n ^2(^i, ^2)^ and notice 
that ri and r2 are treated symmetrically in Eq. (p2D. In what proceeds we will not in detail 
evaluate the integral of every term in this expression. That is a tedious job, and not very 
informative. Instead, we will show an example that is representative and hence give the 
idea for how to do the other calculations. So let us consider the first term in the second 
line of Eq. (P^), and integrate it the way we need in order to find FLe(-R). This yields 



4:71^ driTi / dr2r2{rl - 1) — — [rlg2) 
Jo Jo ri ori ^ ' 

r 1^'!=^? 
47rM dr2r2(r^-l) r^fi'2(ri,r2) 
Jo L 



4vr^^^^/^''dr2r2/.2(i?,r2) 
2 R fdFiR^,R2)^ 



The last equality sign shows the relation between h and one of the derivatives of F. This 
expression easily follows from Eq. (|90|). Other derivatives can be similarly expressed, and 
we use the relations 

= 47.'RiR2kRu R2), (94) 

0K1K2 
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aRi Ri dRi ^ 'A " ' aRi • ' ' 

along with trivial extensions to ^ gj^^/if ^'^WW^' Performing the integrations of 
all terms in Eq. (|92D and collecting together similar terms we end up with the expression 
in Eq. (|3^), that is 



FUR) = FiR,R) + ciiR) 



OF dF 



dRi dR2 



R 



[9^^ d^F\ f d^F 

+C2(^)(l^ + ^) +cl{R) 



R 



dRi 9i?2 / ^ \dRiR2^ 
d'^F d^F \ 2 
dR^^^dKMjR^ \dRlRlj ^ 



+UR)c2iR) ( + +cliR)(,^] . (96) 
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